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Dirac equation for membranes 

M. Trzetrzelewski'Q 
(Dated: January 14, 2013) 

Dirac's idea of taking the square root of constraints is applied to the case of extended objects 
concentrating on membranes in D — 4 space-time dimensions. The resulting equation is Lorentz 
invariant and predicts an infinite hierarchy of positive and negative masses (tension). There are no 
tachyonic solutions. 



MOTIVATION 

In the Dirac membrane model [T] the hamiltonian cor- 
responding to the spherically symmetric solution includes 
a square root of the anharmonic oscillator. As noted by 
Dirac to introduce the spin into the theory one would 
have to get rid of the square root by bringing in the 
spin matrices. However, in the vast literature on mem- 
branes (see e.g. [2] and references therein) a different 
way of introducing fermions for extended objects is usu- 
ally considered by writing the supersymmetric version 
of the bosonic action [3] . Besides aesthetic reasons this 
seems quite not in the direction indicated by Dirac. It is 
perhaps worth mentioning that for the case of point par- 
ticles the Klein-Gordon equation has its supersymmetric 
counterpart (i.e. the massive Wess-Zumino model) which 
however does not describe the electron unlike the square 
root of the Klein-Gordon equation. 

In this paper we choose a route indicated by Dirac i.e. 
introduce fermion fields (and hence spin) for membranes 
by performing a certain square root. We find it more ap- 
propriate to do so for the constraints of the theory rather 
then for the hamiltonian since the former correspond to 
the mass shell constraint for point-like particles and are 
written in relativistically invariant way. However, taking 
the square root of the hamiltonian is also possible and 
was discussed in dj. 

The model presented here is manifestly Lorentz invari- 
ant hence there are no anomalies originating form spe- 
cial relativity. There are also no tachyonic solutions - we 
show that in the case of a spherical membrane the spec- 
trum is real and discrete. The eigenvalues are identified 
with the membrane masses and can be positive as well 
as negative corresponding to positive respectively neg- 
ative membrane tension. The field content (and their 
equations) considered here is novel: the membrane is de- 
scribed by a fermionic field governed by the analog of 
the Dirac equation. It is an alternative way of describing 
extended objects in space-time. 

Lastly we generalize the discussion by including a cou- 
pling to a three form as well as suggest an alternative 



way of writing the Dirac equation for membranes. 



A SQUARE ROOT 

A Lorentz invariant action describing membranes is 
given by Q] 

Smembrane = ~ A / VGda 3 , G a p := daX^dfjX^ (1) 

where A is the tension (here we assume A > and will 
comment on A < later), a a — (r, <r r ) is the internal 
parametrization of the membrane, are the embed- 
ding variables and G is the determinant of the induced 
metric G a p. The momenta := dC/dX^ satisfy the 
constraints (see e.g. [3]) 

p^V" = A 2 det G rs (2) 

and follow from the diffeomorphism invariance of 0. In 
addition to |2]) we also have V^drX^ = 0. 

Equation ^ is a counterpart of the mass shell con- 
straint for a point-like particle therefore in order to keep 
the analogy with the way in which Dirac obtained equa- 
tion for fermions it is desirable to perform the square root 
of(2j. 

Spherically symmetric membrane 

Let us now consider a concrete example of a sphere at 
rest 

X a = t, X 1 = r(r) sin y> cos 0, X 2 = r(r) sin ip sin 0, 

X 3 = r(r) cos<y9, tpe[0,n], 6>e[0,27r). 

The determinant in ^ is det G rs — r 4 sin 2 <p hence the 
corresponding square root is 

7^ = -AAr 2 sin <p (3) 
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where A is a matrix s.t. A 2 — 1, 7^ are gamma matrices, 
{7^,7^} = 2rj^l, the minus sign is a convention. 

Equation pj) is (f dependent. A way to obtain the 
analog of the Dirac equation for the radial variable is to 
average ([3| over the membrane surface by introducing 
the canonical momenta ir u for the whole membrane via 



'ij, ■ — 1 VpdtpdO. 

J T— const 

This implies 

7 m tt m = —M, M := Ar 2 / A sin tpd(pd6 

and suggests to consider the Dirac equation (substituting 
7T„ = —id^) of the form 



(4) 



where ip = ip( r )- 

Because the spinor ip is independent of 9 and ip we may 
consider only the radial part of the differential operator 
in Q. The matrix A should be chosen such that M is 
proportional to the unit matrix. For A = 1 (other choices 
of A are possible which after averaging lead to M = 0), 
going to the radial equation (we use the conventions of 
[6]) we obtain 

(d x + -) G = (e + x 2 )F, f-d x + -)F=(e- x 2 )G 

(5) 

where x — r(47rA) 1 / 3 and e = _E/(47rA) 1 / 3 are dimen- 
sionless, E is the energy, k = ±1 and F(x), G(x) are the 
radial parts of the spinor ip [5]. Using matrix notation 
([5]) can be written as 



Hcf> = ecf>, H := 



The spectrum of H is independent of the choice of n 
(i.e. if <p T = (F, G) solves Hip = e0 for k = 1 then 
<p T = (G,F) solves H(j> — e(f> for k = —1) hence we 
choose k — 1 from now on. Let us note that if <j) is 
the eigenvector of H then (due to the 1/x term) F(0) = 
G(0) = 0. Therefore we may consider the Hilbert space of 
square-integrable vectors 4> on [0, 00) satisfying </>(0) = 0. 
In such Hilbert space the operator d x is antihermitian 
hence H is hermitian which proves that the spectrum of 
H is real. 

Let us also observe that H 2 can be written as 





2/x 2 



Q 



-x 2 d x 
-dx x 2 



hence the spectrum of H 2 is discrete, due to the inequal- 
ity H 2 > Q 2 and the fact that Q 2 is discrete (to see 
explicitly that Q 2 is discrete it is useful to introduce 



a± := F ± G for which the eigen-equation Q 2 (j> = rf<f) 
gives h±a± = rj 2 a±, h± :— —d 2 + ± 2x, i.e. a+ and 
a_ decouple and since /i+and /i_ are discrete, Q 2 must 
also be discrete). This proves that the spectrum of H 2 
and hence H is discrete. 

To find the exact spectrum of H we use numerical 
methods (see the Appendix for details). The first pos- 
itive energy levels are 

ew2.7, 4.0, 5.1, 6.1,7.0, 7.9, 8.7, 9.5, ... 

while the negative ones are 

-ewl.7, 3.3, 4.5, 5.6, 6.6, 7.5, 8.3, 9.1, .... 

The four degrees of freedom of the spinor ip should 
be interpreted just like in the Dirac equation i.e. par- 
ticles and antiparticles with spin (the notion of spin is 
not related here to any geometrical configuration such 
as e.g. the spinning membrane). The antiparticle cor- 
responds to the negative energy — E however the energy 
in this membrane model is a mass. Therefore one con- 
cludes that the model predicts the existence of particles 
with negative masses. Because the energy is measured in 
units of A 1 / 3 and since (—A) 1 / 3 = —A 1 / 3 we can iden- 
tify the negative masses with the negative tension - a 
possibility already considered by Dirac (see the second 
reference in £Q). Such objects are unstable (unless one 
considers gravitational effects with extra volume term in 
the action [T]). The inclusion of the electromagnetic field 
would only worsen the stability of particles with negative 
tension. 

The electromagnetic field can be added using minimal 



coupling prescription Vu —> Vn + eA^ with A 



-e/2r 



(the electrostatic energy of the Coulomb field of a charged 
sphere) and Ai — 0. Equations ^ stay the same with 
shifted energy e — > e + ^- where a is the fine structure 
constant. The corresponding hamiltonian is 



H 



EM 



H 



2x~ 



Due to the small value of a the spectrum of H is 
very close to the spectrum of H (differences are of order 
10~ 3 , the lowest positive energy is now 2.7570 compared 
to 2.7525 obtained previously). 

Let us compare the spectrum of ti EM with the spec- 
trum of the hamiltonian for the bosonic membrane ob- 
tained by Dirac pQ 



H Dirac = V-d? + + 



2r : 



(6) 



where u> — Air A. Using the Bohr-Sommcrfeld quantiza- 
tion method, Dirac found that the first excitation of the 
membrane corresponds to the energy w 53m e , where m e 
is the mass of the electron. The exact value turns out 
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to be about 43m e as observed in [7]. Using the relation 
u = 16ml /27a 2 [1] we find that our model gives 



E= (4 7 rA) 1 / 3 e 



1 f4 



3 V a 



2/3 



m e • 2.757 = 61.5m e 



Therefore we obtained a larger value - in accordance with 
Dirac's expectation that bringing spin into the theory 
might lift the eigenvalue. However the result is still not 
of order of the muon mass. 



General case 

Taking the square root of ^ is more difficult. While 
the l.h.s. would simply give J^V^ the r.h.s. involves a 
square root of the determinant det G rs which is 4th order 
in derivatives. It turns out that the answer is simple in 
terms of a hermitian mass matrix density 



M := -iAe rs 7fiU d r X^d s X v 



(7) 



(where 7 M!/ = Y^Vp^cru, 1^ ■= e rs is com- 

pletely antisymmetric with e 12 = 1) due to the identity 



M 2 = A 2 det G rs l 



(8) 



which allows to write the square root of the constraint as 



= -m. 



(9) 



It is possible to introduce a matrix A such that the rescal- 
ing Ma '■= A - M. is consistent with the constraint 
One simple example is A = AM. / AV det G rs where ma- 
trix A as s.t. A 2 — 1, for which Ma = AA^detG rs l. 
Therefore this particular choice of A will reproduce the 
results of previous section. 

The choice A = 1 results in the averaged mass matrix 
M equal for spherical membrane, hence the spectrum of 
the corresponding Dirac operator would be continuous. 
We now elaborate on this case in more details since it 
is where quantization procedure is possible in terms of 
matrix regularization [5]. 

If the parameters a r are not integrated out the co- 
ordinates X 1 * are fields and the differential equation of 
the quantum counterpart of ([£]) is a functional one. Be- 
cause Vn and X M are conjugate variables it is natural 
to substitute with a functional derivative — ijj^i in 
|9]) . A more problematic term is the mass matrix density 
M. Observing the appearance of the Poisson bracket 
{X^, X v } = e rs d r X»d s X u in Q we write 



5X» 



-iA 7flu {X^,X u })^ = 



(10) 



which will be useful when using the matrix regularization. 
The above equation is similar to Eq. (47) of Ref. [9] 
where an attempt to quantize the bosonic membrane in 
a covariant fashion was undertaken. 



Equation ( 10 ) at least formally is solved by 



e 5 *o, S=^ up j X»{X»,X<>}d 3 a (11) 



where ^0 is a constant spinor, j^p — 5 (7^7^ + cycl.) 
and where we used j^j^p = 2 7l/p . Similar functional 
expressions were found for the bosonic membrane by 
Smolin 9 , Moncrief [TU] and Hoppe [TT]. Using the 
identities 7^ = -707^70 and e s = 7oe _S 7o it fol- 
lows that — ^o^o therefore W suffers from being 
non-normalizable (i.e. the functional integral J[dX}^> , & 
is infinite). 

The hamiltonian formulation can be obtained by sin- 
gling out the hamiltonian operator % := Vq and noting 
that the gauge choice X° — cr° for the coordinates implies 
that {X , X v } = hence 

m = Ef, H = 7 [^Vi + Iay 3 { x \ x3 }) ■ (12) 



In addition to ( 10 1 we also have the constraint "P M <9 r AT A 



which after setting the gauge X L 



a becomes 



VidrX' 1 = 0. Locally (on the membrane surface) this 
constraint is equivalent to {Vi,X 1 } = which will be 
further used. One can write a solution analogous to ( 11 ) 
corresponding to E = 



* = e s V , S = ^ 7 y fc j X l {Xi,X k }d 2 a 



(13) 



(where we used Ylijk — Ijk) which however is not nor- 
malizable by similar argument as before. 



Using matrix regularization [8] we can view ( 12 ) as the 
N — > 00 limit of 

H«f = EM, H = 7° 7 > i + i A 7 ° 7lJ [x\ x 1 } 

where pi and x l are su(N) matrices, subject to the con- 
straint [pi,x l ]^> = 0. The solution corresponding to 
E = is 



V = e b V , S 



A 



6V 7 N 2 - 1 



ll]k Tr{x l [x\x k ]) (14) 



(using the conventions of [H]). The results of the previ- 
ous section suggest that the spectrum of membrane exci- 
tations does not contain the state with E = therefore 



one should not worry that the wavef unctions ( 13 1 and 



( 14 ) are not normalizable. However the exponential fac- 



tor e s will probably play an important role in search for 
normalizable excitations. In the matrix-regularized case 
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this would imply that the tension has to be renormalized 
as A - TV. 

The electromagnetic field can be introduced using 
the minimal coupling prescription — ij^r — > ^ *3xjt + 
eA p . How one can write the matrix regularized case is 
less obvious since A p are in general non polynomial. 



3-FORM COUPLING 

Let us consider a different approach based on p-form 
electrodynamics |13j which is more natural for extended 
objects withp— 1 spatial directions. We therefore assume 
the existence of a three form A in space-time, with com- 
ponents given by a completely antisymmetric tensor Ap, vp 
and the corresponding field strength T pvpa = d^A^pa] 

For the action of the membrane coupled with A we take 
S = - A J ^Gd 3 a-^ J t ap ''A^ p d a X p d p X v d 1 X p d 3 a 

so that the momenta and the equations of motion are 
V p = - A ^/\G\G T|3 dp X p - ^eA„ vp {X»,X>>}, 

Ad a VGG a0 d p X p + eT p vp<y {X\X p , X a } = 

where {X p ,X V , X p } := e a ^ d a X p dpX v d^Xf is the 
Nambu bracket [14] . The corresponding square root of 



constraints will result in (|10|) with the functional deriva- 
tive -ij^ -> ■ " 



'6» 



^eA^plX" , X p }. One can also 
write the matrix regularization. 

Let us note that the equations of motion for A pvp do 
not contradict sourceless Maxwell equations in the fol- 
lowing sense: using the identity F^vpo — —6(- P v P trd 7r a 7r 
where a lr :— \t K>11 ' p A llV p {t pup<y is completely antisym- 
metric, £i234 = 1), the action for J- becomes 



S 



1 



2 • 4! 



T^ pa T plJpa d A x = 18 / {d^fd^x 



i.e. a square of the Lorenz gauge condition. The equa- 
tions of motion d v {d p a p ) = imply d v a v = const, hence 
the corresponding space of solutions is larger then that 
of the electromagnetic potential A p in the Lorenz gauge. 
Moreover the residual gauge invariance of T 2 is present 



with a,. 



d P Xi — just like in the formulation 



of electrodynamics in the Lorenz gauge. 



ALTERNATIVE SQUARE ROOT 

An identity similar to ^ holds also for the whole 
world- volume metric i.e. 



Both identities have also a heuristic origin when intro- 
ducing a matrix line element dX := j^dx^ and noting 
that components of the matrix surface and volume ele- 
ments are 

dX AdX = ^{X^^X^jda 1 A da 2 

for one forms dx p — d r X p da r on the surface and 

dX AdX AdX = ^^{X p , X",X p }da Q A da 1 A da 2 

for one forms dx^ — d a X p da a on the world-volume. 
Therefore matrices Ai and C appear quite naturally when 
working with matrix line element dX. 

Having in mind the above remarks let us consider the 
following analogy. In the case of point particles the La- 
grangian formally can be written as a matrix 

yJx^dT = j^dr = ■y p dx p 

which has a counterpart in —i"/ p d p - the differential part 
of the Dirac operator. Likewise, for membranes we have 

VGd 3 a = l - lp „ p {X»,X»,X p }d 3 <r = i lpvp dn pvp 

where dfl pup := \{X P , X v , X p }d 3 a is the element of the 
world-volume, hence the corresponding Dirac operator 
would be "f pi/p fiQ^p , with the functional derivative (al- 
ready considered by Nambu |13j). Therefore in the pres- 
ence of the field A pvp we would obtain 



r 



-+eA^ p U> = 



which is an alternative way to write the Dirac equation 
with the three-form. 



CONCLUSIONS 

A theory of extended objects aspiring to describe par- 
ticles should include fermionic fields. How to write down 
such theory is non-trivial considering the fact that the 
starting point is a classical object described by bosonic 
fields (coordinates) . A problem of this kind but for point- 
like particles was solved by Dirac by taking a square root 
of the constraints and since it is the correct route we 
take it as a guiding principle in the attempt to include 
fermions into the theory of extended objects. 

Because the constraints are local equations depending 
on the geometry of the membrane, the corresponding 
Dirac equation will be a functional one. Even for the 
simplest choice of the matrix A, this results in a highly 



C z = detG a/3 l 



C 



3! 



1 ^ p {X p ,X\X p }. 



complicated system ( 10 1 which fortunately can be under- 



stood using the matrix regularization so that the masses 
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of the membrane are hidden in the spectral properties of FIG. 1: Convergence of the energy levels with the cutoff up 
such matrix-Dirac operator. Still to find the eigenvalues to n — 300. 
of this equation is probably very difficult which is why 
we considered a semiclassical approach by concentrating 
on spherically symmetric membranes. The resulting sys- 
tem is a Dirac equation for the radius of the membrane. 
Although this approach is only approximate we believe 
the we have identified the crucial properties of the theory 
i.e. that the spectrum is discrete and consists of positive 
as well as negative eigenvalues - in particular there is 
no eigenvalue and there are no tachyons. The negative 
masses are due to the negative tension and therefore such 
particles would be unstable. 
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Appendix 

In order to obtain numerically the eigenvalues of matrix 
operator H we use the cutoff method |15) which consists 
of calculating the representation of H in some orthonor- 
mal basis (we choose , 1 ^=xL^ (x)e~ x ^ 2 , where 

V V(n + l)(n+2) W 

L^n\x) are generalized Laguerre polynomials), truncate 
the infinite matrix (i.e. truncate each operator appearing 
in the entries of the 2x2 matrix H) and then numerically 
diagonalize it. The spectra and the eigenvectors of the 
truncated matrices converge to their exact counterparts. 
The results of the numerical approach are presented in 
Figure 1. 
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